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After a brief historical introduction to Bose-Einstein condensation and Fermi de- 
generacy, we discuss theoretical results we have recentely obtained for trapped 
degenerate quantum gases by means of a thermal field theory approach. In par- 
ticular, by using Gross-Pitaevskii and Bogoliubov-Popov equations, we consider 
thermodynamical properties of two Bosonic systems: a gas of Lithium atoms and 
a gas of Hydrogen atoms. Finally, we investigate finitc-tempcrature density profiles 
of a dilute Fermi gas of Potassium atoms confined in a harmonic potential. 



1 Introduction 

Bose-Einstein condensation (BEC) is the macroscopic occupation of the lowest 
single-particle state of a system of bosons. The concept of Bose condensation 
dates back to 1924. In that year it appeared a paper of Bose [1] about a 
new derivation of photon statistic and Planck distribution and the same year 
Einstein [2] gave the theoretical description of BEC for a homogeneous system 
of identical particles. Einstein used the de Broglie matter wave for applying 
Bose calculation and developing a statistic also for particles. In 1938 London 
proposed the idea of a macroscopic wavefunction to describe the behavior 
of a Bose-condensed fluid [3]. His idea was that BEC should be involved in 
the phase transition of superfluid He. From this suggestion was derived the 
Tisza's two- fluid concept [4] and, consequently, the Landau's model [5]. Other 
contributions came from the Russian school (Bogoliubov, Beliaev, Popov). 
In particular the study of BEC and its elementary excitations in a weakly 
interacting dilute Bose gas [6] . 

From 1995 we have strong experimental evidence of BEC in clouds of 
confined alkali- metal atoms at ultra-low temperature (about 100 nK): at JILA 
with 87 Rb atoms [7], at MIT with 23 Na atoms [8], and at Rice University 
with Li atoms [9] . In these experiments the condensed fraction can be more 
than 90%. More recently, another MIT group has obtained BEC (but small 
condensed fraction) with a gas of hydrogen atoms [10]. 

For dilute vapors one can separate off the effect of BEC and study the 
role of the interaction. The theoretical understanding of recent experiments is 
based on the Gross-Pitaevskii equation of a space-time dependent macroscopic 



wavefunction (order parameter) that describes the Bose condensate [11]. 

The experiments with alkali-metal atoms generally consist of three steps: 
a laser cooling and confinement in an external potential (a magnetic or 
magneto-optical trap), an evaporative cooling and finally the analysis of the 
state of the system. It is important to observe that the Nobel Prize 1997 has 
been given for the development of methods to cool and trap atoms with laser 
light [12]. Nowadays more than twenty experimental groups have achieved 
BEC by using different geometries of the confining trap and atomic species. 

It is interesting to note that the paper of Einstein on BEC [2] preceded 
the concept of Fermi-Dirac statistics [14] as well as the division of all particles 
into two classes (Bosons and Fermions) depending on their net spin [15]. 

Very recently the Fermi quantum degeneracy has been obtained with di- 
lute vapors of laser-cooled 40 K atoms confined in a magnetic trap [15]. This 
experiment has renewed the interest on Fermi gases. In fact, with an effective 
attractive interatomic interaction, it should be possible to produce Cooper 
pairs of Fermions with Bosonic statistics and to investigate the mechanism of 
superfluidity with Fermions. 

In this paper we investigate dilute quantum gases at ultra-low tempera- 
ture and confined by an external potential. In section 2 we review some basic 
concepts of thermal field theory for Bosons. In sections 3 and 4 we apply 
thermal field theory to study a Bose gas of Lithium and Hydrogen atoms, 
respectively. In section 5 we discuss thermal field theory for Fermions. Fi- 
nally, in section 6 we analyze some properties a dilute Fermi gas of Potassium 
atoms. 



2 Thermal Field Theory for Bosons 

Thermal field theory investigates the many-body dynamics of finite- 
temperature systems, whose underlying structure is represented by quantum 
fields. This theory is the combination of statistical mechanics and quantum 
field theory [16]. Here, we apply the thermal field theory for non-relativistic 
Bosons in the same hyperfine state. 

The Heisenberg equation of motion of the bosonic field operator ip(r,t), 
which describes a non-relativistic system of confined and interacting identical 
atoms in the same hyperfine state, is given by 



ih-j(r,t) = [-— V 2 + [/(r)]^(r,i) + J rf 3 rV+(r',i)y(r,r')^(r',^(r,t) , 

(1) 

where m is the mass of the atom, U(r) is the confining external potential and 
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V(r, r') is the interatomic potential. 

The Bosonic field operator must satisfy the following equal-time commu- 
tation rules 



(2) 
(3) 



[^ + (r',t),^(r,t)]=^(r'-r), 

mr',t)^(r, t)} = $+(r>, t),j+(r, t)]=0, 

where S 3 (r) is the 3D Dirac delta function and [A, B] = AB — BA. 

When the density of the atomic cloud is such that the scattering length 
and the range of the interatomic interaction are less than the average inter- 
atomic distance, the true interatomic potential can be approximated by a 
local pseudo-potential 

V(r,r')=.g<5 3 (r-r'), (4) 

where g — 4irh 2 a/m is the scattering amplitude of the spin triplet channel (a 
is the s-wave scattering length). In this way the Hamiltonian operator H of 
the system reads 

H= fd 3 r Vi+(r) [ - |-V 2 + U(r)] $(r) + ^W> + (r)^ + (r)^(r)^(r) . (5) 



Because the Hamiltonian H is invariant under the global U(l) field transfor- 
mation 

$(t) e">(r) , (6) 
the Noether theorem gives us the conservation of the particle number 

N = J d 3 r V»+(r)^(r) . (7) 

In equilibrium satistical mechanics, the grand canonical partition function 
Z of a quantum system of Hamiltonian H and conserved number N of particles 
is written as 

~ e -{H-iiN)/kT~ 



Tr 



(8) 



where T is the temperature of the thermal reservoir, k is the Boltzmann 
constant, and \i is the chemical potential [16]. The thermal average of an 
observable A is given by 



(A) = -Tr 



A e -(H-fiN)/kT 



(9) 



Note that the chemical potential fj, fixes the average total number of parti- 
cles. In the grand canonical ensemble the system is described by an effective 



Hamiltonian H — fiN. Therefore, the chemical potential n can be introduced 
with the following shift 

d d i ,„„, 

in the equation of motion of the field operator ^(r, t). In this way the equation 
of motion of the field operator at finite temperature is given by 

► 2 



d 



2m 



iK-ilj{r,t)= -—X7 2 + U(r)-LiiP(r 7 t)+g^+(r,tW(r,t)^{r,t). (11) 



In a Bosonic system one can separate Bose-condensed particles from non- 
condensed ones by using of the following Bogoliubov prescription 

$(r,t) = *(r)I + 0(r,t) (12) 

where / is the identity operator and 

*(r) = $(r,t)) (13) 

is the order parameter (macroscopic wavefunction) of the condensate. The 
field <l>(r,t) is the fluctuation operator, which describes the non-condensed 
fraction of Bosonic atoms [6] . It is important to observe that the Bogoliubov 
prescription breaks the U(l) symmetry of the system, namely the number 
operator N is not a conserved quantity if <f>(r) ^ 0. For a fixed chemical 
potential /u, it exist a critical temperature Tbec, called BEC critical temper- 
ature, below which the order parameter is no more identically zero [17] . In this 
case one speaks of spontaneous symmetry breaking because the ground-state 
does not possess the initial symmetry of the system [18]. This phenomenon 
is also known as off diagonal long-range order because, in presence of BEC, 
the two-body density operator (^(r')^(r)) must satisfy the Penrose- Onsager 
criterion 

lim (VS(r')VKr)) = $*(r')*(r) , (14) 

|r'— r| — >oc 

where $(r) is the eigenfunction of the largest eigenvalue of the two-body 
density operator [19]. 

The Bogoliubov prescription for the field operator ip(r,t) enables us to 
write the three-body thermal average in the following way 

(^+(r, t)tjj(r, t)tjj(r, t)) = |$(r)| 2 $(r) + 2n(r)$(r) + m(r)$*(r) + s(r) , (15) 

where n(r) = (0 + (r, i)</>(r, t)) is the density of non-condensed particles, while 
m(r) = (cj)(r,t)(j)(r,t)} and s(r) = (0 + (r, t)(f>(r, i)</>(r, t)) are anomalous den- 
sities. Now, we obtain an equation for $(r) by taking the thermal average in 
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the equation of motion for the field operator 4>(r). In this way we find 



2m 



V 2 + J7(r)$(r)+#|$(r)| 2 + 2 5 n(r) $(r)+ ff m(r)$* (r) + S(r) = /z$(r) , 

(16) 

that is the exact equation of motion of the Bosc-Einstein order parameter 
<5>(r). This is not a closed equation due to the presence of the non-condensed 
density n(r) and of the anomalous densities m(r) and s(r). Obviously, if the 
particles of the system are all in the Bose condensate then the non-condensed 
density and the anomalous densities are zero. In this case the previous equa- 
tion becomes 



2m 



V 2 + C/(r)$(r)+ ff |$(r) 



/U$(r) 



(17) 



that is the so-called Gross-Pitaevskii equation [11]. This equation can be used 
to describe a dilute Bose gas near zero temperature, when the non-condensate 
is very small compared with the condensate fraction. Another possible, less 
drastic, simplification, neglects only the anomalous densities. Then the equa- 
tion of motion of the Bosc-Einstein order parameter $(r) becomes 



2m 



V 2 + f/(r)$(r) + ,g|$(r)| 2 + 2gn(r) $(r) = /z$(r) . (18) 



Also this equation, that we call finite-temperature Gross-Pitaevskii equation, 
is not closed. We must add an equation for the non-condensed density n(r) 
by studying the fluctuation operator (j)(r,t). 

The exact equation of motion of the fluctuation operator 0(r, t) is ob- 
tained by subtracting Eq. (16) to Eq. (11) and using Eq. (12). Such equation 
is given by 



|-V 2 + ^) 



M + 2.g|$(r)| 2 0(r,t)+ 



+ g <D(r) 2 0+(r,t) + <f>*(r)0(r,t) 2 + 2$(r)^+(r,i)0(r,i)+ (19) 



+0+ (r, f)$(r) 2 0(r, tf - m(r)$(r) - s(r)$(r) 

To find the thermal excitations above the condensate one can simply linearize 
the previous equation. This can be done by means of the mean-field approxi- 
mation: 



4>+(T,t)4>(T,t) ~ (4>+(T,t)4>(T,t))=n(T) 



(20) 
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^(r,t)^(r,t) ~ {(j>{Y,t)cj){v,t)) = m(r) , 



(21) 



0+(r,t)0(r,t)^(r,i) ~ 2(<j>+ (r,t)<j>(r,t)}<j>(v,t) + (0+(r,t)^(r,t)^(r,t)) = 



= 2n(r)0(r, t) + m(r)0+ (r, t) + s(r) 



(22) 



Within the mean-field approximation, the linearized equation of motion of the 
fluctuation operator reads 



— V 2 + U(v) M + 2. 9 (|<I>(r)| 2 + n(r))p(r, i)+ 



g($(r) 2 + m(r))^ + (r,t) 



(23) 



This equation remains very complex and some simplification is usually per- 
formed. One possible approximation, called Bogoliubov standard approxima- 
tion, neglects the non-condensate densities (n(r) = m(r) = 0). This ap- 
proximation is particularly useful near zero temperature: the condensate and 
its elementary excitations are separately calculated. Another possible, less 
drastic, simplification is called Popov approximation. It neglects only the cor- 
relation function m(r) and it is able to give a more accurate description of the 
non-condensate fraction. The Popov approximation is expected to be reliable 
in the whole range of temperature except near T c , where mean-field theories 
are known to fail. 

In any case, due to the presence of the self-adjoint operator cf) + (r, t) in the 
linearized equation of motion (23), the fluctuation operator must be expanded 
by using the so-called Bogoliubov transformation 



4>(r,t) = J2 [u j (r)e- iE ^ h a j +v*(r)e iE ^ h a+ 



(24) 



where &j and are bosonic operators and the complex functions Uj(r) and 
Vj (r) are the wavefunctions of the so-called quasi-particle and quasi- hole exci- 
tations of energy Ej . The functions Uj (r) and Vj (r) satisfy the normalization 
condition 



/ 



d 3 r [u*(r)u k (r) - v*(r)v k (r)] = 5 jk . 



(25) 



By inserting Eq. (24) into Eq. (23) with m(r) = 0, we get the Bogoliubov- 
Popov equations 

h 2 

~ 2^ V2 + U{r) ~ M + 2gn{r) \ %(r) + 5<f>(r)2 ^ (r) = E ^ r) ' 
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■ £ y2 + u[r) ~ fi + 2gn[r) ] v ' j[v) + 5<f> w 2 ^ r ) = - e ^ t ) < ( 2g ) 

where n(r) = |$(r)| 2 + n(r) is the total density. 

Eq. (18) and Eq. (26) are supplemented by the relation fixing the total 
average number of atoms in the system 

N = J d 3 r[n (r)+n(r)] , (27) 

where 

n (r) = |<f>(r)| 2 , 
n(r) = (l^«| 2 + l^«| 2 ) + I^WI 2 > (28) 



with 



1 



W = e E j/ kT_ 1 ( 29 ) 

the Bosc factor at temperature T. Note that also at zero temperature there 
is a non-condensed density (quantum depletion), given by J^. |fj(r)| 2 . 

When fcT is much larger than the lowest elementary excitation, one can 
use the Hartree approximation neglecting the quasi-hole excitations (vj (r) = 
0) . In this way, the Bogolibov-Popov equations reduce to the following Hartree 
equations 

K 2 1 
— V 2 + U(r) - n+ 2gn(r)\uj(r) = E jUj (r) , 

and the non-condensed density is simply 

3 

In this case the non-condensed density is equivalent to the thermal density 
because the quantum depletion is zero. If there is a large number of particles 
in the thermal cloud, the Hartree thermal density can be calculated through 
the quasi-classical approximation. In the quasi-classical approximation one 
uses the classical single-particle energy E(p,r) = p 2 /(2m) + U(r) + 2gn(r) 
instead of the quantum energy Ej . Thus the thermal particles behave as non- 
interacting Bosons moving in the self-consistent effective potential U(r) + 
2gn(r). Then the previous formula becomes 

d3 P * = J_ „„,„ /V(E/(r)+2 S n(r)-^)/fcT^ 

(31) 



, f d^p 1 1 / 

nW = J (2^)3 e(^)-ri/^-l = 53/2 I" 
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where A = (2ttTi 2 /mfcsT) 1 / 2 is the thermal length and 



oo 



k 



9 a {z) = 



z 



(32) 



k=i 



k a ' 



In the effective potential, the term 2gn(r) is the Hartree mean-field generated 
by interactions with other atoms. 

3 Bose gas of Lithium atoms 

In this section we study the thermodynamics of a Bose-Einstein condensate 
of 7 Li atoms, which is particularly interesting because of the attractive inter- 
atomic interaction: only a maximum number of atoms can form a condensate; 
beyond that the system collapses. Up to now there are only few experimental 
data supporting BEC for 7 Li vapors with a limited condensate number. Be- 
cause of the small number of condensed atoms, the size of the Bose condensate 
is of the same order of the resolution of the optical images. So it is difficult to 
obtain precise quantitative informations, like the number of condensed atoms 
and the BEC transition temperature [9]. 
We consider an isotropic harmonic trap 



and use trap parameters of the experiment [9] : lu = 878 Hz for the frequency of 
the isotropic harmonic trap and a = — 27ao for the scattering length (ao is the 
Bohr radius). We measure the radial coordinate r in units of the characteristic 
oscillator length ah = \/%/ (mui) and the energies in units of energy oscillator 
quantum fiuj. 

We concentrate our attention on density profiles and condensate fraction 
by using both the Bogoliubov and Popov approximations. Actually, the Popov 
method and its quasi-classical approximation are a good starting point to 
estimate the contribution of the correlations in the gas for any finite value of 
the temperature. 

In the case of negative scattering length, a simple variational calculation 
with a Gaussian trial wavefunction shows that a harmonic trap supports a 
Bose condensate with a number of bosons smaller than a iV c = 0.67(o/j/|a|) 
[21]. In a homogeneous gas such critical number is zero. 

The system with TV = 1000 atoms is sufficiently far from the critical 
threshold N c = 1260. We solve self-consistently the finite-temperature Gross- 
Pitaevskii equation (18) and the Bogoliubov-Popov equations (26) for some 




(33) 
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Figure 1. Popov calculations for 10 3 7 Li atoms. Profile of the condensate density for tem- 
perature ranging from to 58nK. 



value of the temperature. In this way, we obtain informations about the den- 
sity profiles of the condensate and non condensate fractions and the spectrum 
of elementary excitations [20]. 

We plot in Fig. 1 the profile of the condensate fraction for some value 
of temperature from to 58nK. Note that the temperature does not modify 
the shape of the condensate density and its radial range. This is not the case 
of the non condensed density profile, for which the increasing temperature 
produces a sharper maximum at a larger value of r. 

In Fig. 2 we plot the condensate fraction N /N versus the temperature T. 
In this figure we include also the quasi-classical approximation of the Popov 
and Hartree method and the curve of the ideal gas in the thermodynamic 
limit. We observe that the quasi-classical Popov approximation, for which the 

classical energy reads E(p, r) = \J [p 2 / (2m) + U(r) + 2gn(r)] 2 + g 2 n (r) (see 
also [20]), gives accurate results for the condensate fraction at all temperatures 
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Figure 2. Condensate fraction vs temperature for 10 3 7 Li atoms. 



here considered. 

Now we analyze a total number of atoms near the critical threshold 
N c ~ 1260, at which there is the collapse of the condensate. We use the 
quasi-classical Popov approximation to evaluate the critical threshold iV c as 
a function of the temperature. 

In Fig. 3 we plot the condensate fraction as a function of the temperature 
T for different numbers of 7 Li atoms. By increasing the temperature it is 
possible to put a larger number of atoms in the trap but the metastable 
condensate has a decreasing number of atoms. Also the critical number Nq 
of condensed atoms decreases by increasing the temperature. These results, 
obtained with the quasi-classical Popov approximation, show that, when u) = 
878 Hz, for N < N c — 1260 the system is always metastable and it has a finite 
condensate fraction for < T < T BEC , where Tbec is the BEC transition 
temperature. For N > N c the system is metastable only for temperatures 
that exceed a critical value T c and it has a finite condensate fraction for 
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T C <T<T BEC [20]. 

We have seen that when the condensate fraction is only few percents 
then n(0)/no(0) becomes sufficiently large and the BEC cannot take place 
anymore. Thus, our calculations suggest that there should be a critical tem- 
perature T* beyond which the BEC transition is inhibited, independently of 
the number of atoms in the trap. Obviously, T c and T BEC are functions of N. 
Such temperatures, and also the critical parameter N c , are functions of the 
trap geometry, namely the frequency lu, and of the scattering length a [20]. 

Finally, we observe that our calculations could be useful to select exper- 
imental conditions for which a fraction of condensed 7 Li atoms is detectable 
at thermal equilibrium. 
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4 Bose gas of Hydrogen atoms 

BEC has been also achieved with atomic Hydrogen confined in a Ioffc- 
Pritchard trap [10]. That is an important result because Hydrogen properties, 
like interatomic potentials and spin relaxation rates, are well understood theo- 
retically. The s-wave scattering length of the Hydrogen is very low (a = 0.0648 
nm) and, compared with other atomic species, the condensate density is high, 
even for small condensate fractions. Moreover, due to hydrogen's small mass, 
the BEC transition occurs at higher temperatures than previously observed. 

We analyze the thermodynamical properties of the trapped Hydrogen 
gas by using the quasi-classical Hartree approximation. This approach is 
justified by the very large number of atoms (about 10 10 ) in the trap and by the 
relatively high temperatures involved (order of pK). Our detailed theoretical 
study of the hydrogen thermodynamics can give useful informations for future 
experiments with a better optical resolution and a larger condensate fraction 
[22]. 

In the experiment reported in [10], the axially symmetric magnetic trap 
is modelled by the Ioffe-Pritchard potential 

U(p,z) = v /( a p) 2 + (/3z 2 + 7 ) 2 -7, (34) 

where p and z are cylindrical coordinates and the parameters a, f3 and 7 
can be calculated from the magnetic coil geometry. In particular, for small 
displacements, the radial oscillation frequency is u> p = a/^/my = 2ir x 3.90 
kHz, the axial frequency is uj z = y/2(3/m = 2tt x 10.2 Hz and j/ks — 35 pK 
[10]. 

It is important to observe in the experiments with atomic Hydrogen, the 
gas is dilute (na 3 << 1) but also strongly interacting because Na/au >> 1, 
where ah = (h/mujh) 1 ^ 2 with ojh = (^wj 1 / 3 . It follows that the kinetic 
term of the finite-temperature Gross-Pitaevskii equation (18) can be safely 
neglected and one gets the Thomas-Fermi condensate density 

n (r) = - g [p - U(v) - 2gn{v)] , (35) 

in the region where p > U(r) + 2griT{v) : and no(r) = outside. In practice, 
we study the BEC thermodynamics by solving self-consistently the previous 
equation and Eq. (31). 

In Fig. 4 we compare the BEC transition temperature Tj? of the Ioffe- 
Pritchard potential with the analytic formula T® H = OMTiUhN 1 / 3 /k B , that 
is exact in the thermodynamic limit for the harmonic potential. Fig. 4 shows 
that, for a large number of atoms, T®' H exceeds Tj?. In particular, for TV = 
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i / • loffe-Pritchard potential (interacting) 

- 1 / ■ Harmonic potential (interacting) 

// loffe-Pritchard potential (non interacting) 

f Harmonic potential (non interacting) 



4.0 1 9 8.0 1 9 1.2 1 1 ° 1.6 10 1 ° 2.0 10 1 ° 
N 



Figure 4. BEC transition temperature T c versus number TV of hydrogen atoms. Comparison 
between harmonic and loffe-Pritchard trapping potential. 



2 • 10 10 the relative difference is more than 41%. This strong effect is due to 
the larger density of states in the loffe-Pritchard trap. 

The role of the interatomic interaction on the transition temperature T c 
is very small. In Fig. 4 one observes that the repulsive interaction reduces T c 
both in harmonic and loffe-Pritchard traps. This effect is of the order of 0.1 
fiK. 

In Fig. 5 we plot the condensate fraction N /N as a function of temper- 
ature for TV = 10 8 , 10 9 and 10 10 atoms. As previously stated, such curves 
could be a useful guide for future experiments with a better optical resolution 
and a larger condensate fraction (further details can be found in [22]). 
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Figure 5. Condensate fraction No/N as a function of temperature T. N interacting hydro- 
gen atoms in the Ioffc-Pritchard trap. 



5 Thermal Field Theory for Fermions 

The Heisenberg equation of motion of the fermionic field operator Vv(r,i), 
which describes a non-relativistic system of confined and interacting identical 
atoms in the hyperfine state a, is given by 



3 ; , . 



2m 



V 2 + C/(r) 



+ J2 I ^+(r',t)^,(r,r')i'(r',#,(r,i) , (36) 

a' J 

where m is the mass of the atom, U(r) is the confining external potential and 
V^v(r, r') is the interatomic potential [16]. 

14 



The Bosonic field operator must satisfy the following equal-time anti- 
commutation rules 

[^+(r',i),Vv(r,i)] + = <W 3 (r'-r) , (37) 

[VV (r', t), fa (r, t)}+ = (r', i), ^+ (r, *)]+ = , (38) 

where [A, B] + = AB + BA. 

When the density of the atomic cloud is such that the scattering length 
and the range of the interatomic interaction are less than the average inter- 
atomic distance, the true interatomic potential can be approximated by a 
local pseudo-potential 

V a , a (r,r') = g ala 8 3 (r-r') , (39) 

where g a > a — Airfr 2 a a > a / 'm is the scattering amplitude and a„' a is the s-wave 
scattering length. The s-wave scattering between Fermions in the same hy- 
pcrfine state is inhibited (aw = 0) due to the Pauli principle. Nevertheless, 
the effect of interaction could be very effective for a Fermi vapor with two or 
more hyperfme states. Thus, in the grand canonical ensemble, the equation 
of motion of the field operator, which describes a dilute fermi gas with atoms 
in different hyperfme states, can be written as 

ih—fa{r,t) = [-— V 2 +C/(r)- M ]^(r,i)+^. 9ffV ^+(r,t)^(r,t)^(r,t) , 

(40) 

while in the case of a dilute Fermi gas of atoms in a unique hyperfme state, 
the previous equation is simply 



i) 



h 



2 



jft-VKr, t) = [ - —V 2 + U(r) - M j V(r, t) . (41) 

This is the equation of motion of an ideal fermi gas in an external trapping 
potential. 

For an ideal Fermi gas in external potential the Fermionic field can be 
expandend in the following way 

^(r,i) = ^n J (r)e-^' t /^, (42) 

j 

where Uj(r) is the single-particle eigenfunction with eigenvalue ej and aj is 
the lowering Fermi operator of the single-particle eigenstate The thermal 
spatial density is given by 

n(r) = $+(r,t)^(r,t)) = ^ \u 3 (r)\ 2 (a+a 3 ) , (43) 

3 
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with 

(*J & j) = e E 3 /kT + I ( 44 ) 

the Fermi factor at temperature T. In the quasi-classical approximation one 
uses the classical single-particle energy E(p,r) = p 2 /(2m) + U(r) instead of 
the quantum energy Ej . In this way the previous formula becomes 

n (r) = f -^E - = — f,„ ( e -(U(r)-n)/kT\ (45) 

where A = (2nh 2 /inksT) 1 / 2 is the thermal length and 

1 poo T n-1 

'•"^rwl ^FvTT' < 46) 

with r(x) the factorial function. Note that for \z\ < 1 one has f n (z) = 
YT=i^ l ) l+1 zl I in ■ Moreover, by using g n {z) = -f n (-z) instead of f n (z), 
one finds the spatial distribution of the ideal Bose gas in external potential. 
In the limit of zero temperature, with fj, = Ep the Fermi energy, the quasi- 
classical spatial distribution is given by the Thomas-Fermi approximation 

n(r) - (2m) 3 / 2 /(67r 2 ^ 3 ) {E F - C/(r)) 3/2 6 {E F - U(r)) , (47) 

where & is the Heaviside step function [17]. 



6 Fermi gas of Potassium atoms 

In 1999 the Fermi quantum degeneracy has been obtained with dilute vapors 
of laser-cooled 40 K atoms confined in a magnetic trap [13]. In that experiment, 
to favor the evaporative cooling, a 40 K Fermi vapor in two hyperfine states 
is used. When the system is below the Fermi temperature, one hyperfine 
component is removed and it remains a trapped quasi-ideal degenerate Fermi 
gas. The s-wave scattering between fermions in the same hyperfine state is 
inhibited due the Pauli principle. It follows that at low temperature the dilute 
Fermi gas, in a fixed hyperfine state, is practically ideal. Nevertheless, the 
effect of interaction could be very effective for a Fermi vapor with two or more 
hyperfine states. 

In the case of a harmonic external potential 

U(r) = ^m{cj 2 x 2 + cv 2 y 2 + u\z 2 ) , (48) 
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one finds the Fermi density profile by using the Eq. (43) and the eigenfunctions 
u nin2ri3 (r) of the harmonic oscillator 

n (r) = V \u ni n 2 n 3 (r)\ 2 , . 

W Z^i e /3ft(wi(m+l/2)+a;2(n2+l/2)+a;3(n3+l/2)-f*) _|_ 1 " V > 

Because the Fermi gas is ideal, one has u ni n 2 n 3 (r) = u ni (x)u n2 (y)u n3 (z), 
where <j) n (x) is the eigenfuction of a ID harmonic oscillator with frequency u 
and quantum number n. This eigenfunction can be found by means of the 
recursion relation 

u n (x) = -^= V2a h xu n -i(x) - \fn - lu n - 2 {x) , (50) 

where uq(x) = a 1 l / 2 ir 1 / 4 e~ a2hX ' "'I 2 and u\(x) = y/2ahxu a (x) 1 with ah = 
yjTi/ (mijS). 

In the quasi-lD case, namely a cigar-shaped gas where lo\ = u>2 » ^3, 
the shell effects arc strongly enhanced. In Fig. 6 we plot the density profile 
of a ID ideal Fermi gas in harmonic potential as a function of the number 
N of particles at kT/huj = 1CP 3 . We have verified that at this temperature 
the density profiles coincide with the zero-temperature ones. The results are 
obtained by numerically evaluating expressions (49) and (50) in the ID case. 
The local maxima, whose number grows with N, are clearly visible for a small 
number of particles. 

In Fig. 7 we show the density profiles with 20 and 50 particles as a function 
of temperature. Remarkably, the local peaks are no more distinguishable for 
temperatures well below the Fermi temperature Tp — Nhuj/k. Thus, to 
see spatial shell effect on the ID density profile, the system should be at 
temperatures lower than Tp by one or two orders of magnitude. 

The problem of a dilute Fermi vapor with M hyperfine states (compo- 
nents), can be studied by using the mean-field approximation and the quasi- 
classical formula. The spatial density profile n a (r) of the er-th component in 
a M-component Fermi vapor can be written as 

Mr) = ( e (--^)-£^^^))/^ (51) 

where a = 1, M, [i a is the chemical potential of the c-th component, and 
g a i a = ATTh 2 a a i a /m, with a a i a the s-wave scattering length between a'-th 
and <7-th component (a aa = 0). Thus, the effect of the other M — 1 Fermi 
components on the a-th component is the appearance of a mean-field effective 
potential. 
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We numerically solve the set of equations (51) with a self-consistent it- 
erative procedure. If the components of the Fermi vapor are non-interacting 
then they can occupy the same spatial region. Instead, if they are strongly 
interacting (repulsive interaction) there will be a phase-separation, i.e., the 
Fermi components will stay in different spatial regions. 

When two components have the same number of particles, the onset of 
phase-separation is also an example of spontaneous symmetry breaking. In 
fact, if the chemical potentials /ij of the two components are equal, Eq. (51) 
always admits a symmetric solution ni(r) = n 2 (r). However, for particle 
number N larger than a threshold N c the solution bifurcates and a pair of 
symmetry breaking solutions appears (see [23]). 

Phase-separation also appears in a Fermi vapor with three or more com- 
ponents. In Figure 8 we plot the density profiles of the 40 K Fermi vapor with 
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Figure 7. Density profiles for an ideal Fermi gas in a ID harmonic trap at finite temperature, 
with N = 20 and 50. N is the number of Fermions. a) kT/hu = 1/10; b) kT/hu = 1/2; c) 
kT/hu) = 1; d) kT/huj = 5. Units as in Figure 6. 



three components in a 3D isotropic harmonic potential. In this case we nu- 
merically solve Eq. (51) with an = ai3 = (I23 = a. The Figure shows that 
also for three components the phase-separation and the apparence of spatial 
shells are regulated by the scattering length. 

A Fermi vapor with more components has the same behavior. In par- 
ticular, we have recently shown [23] that the critical density of fermions at 
the origin, which gives rise to the phase-separation, does not depend on the 
number of Fermi component and on the properties of the trap; moreover we 
have found, by using standard bifurcation theory, that such a density is given 
by n c (0) = 7r/(48a 3 ), where a is the s-wave scattering length. 
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Figure 8. Density profiles of the 40 K vapor with three components (solid, dotted and dashed 
lines) in a 3D isotropic harmonic trap (a; = 450 Hz). N = 0.5T0 7 atoms for each component. 
Different values of the scattering length a (ao is the Bohr radius). Temperature T = 0.57V 
(Tp = 1.07 £tK). Lengths in units = (6/mu) 1 / 2 and densities in units <J^ 3 . 



7 Conclusions 



We have shown that non-relativistic quantum field theory is extremely use- 
ful to explain results of recent experiments with dilute alkali-metal atoms 
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at ultra-low temperatures and to predict new quantum phenomena not yet 
experimentally observed. By studying these systems, sophisticated theoreti- 
cal concepts, like order parameter, spontaneous symmetry breaking and off- 
diagonal long-range order, find a clear correspondence in measurable physical 
quantities. Finally it is important to stress the role of the quasi-classical ap- 
proximation, which is not only useful in computations but is also a giude for 
a deeper understanding of physical problems. 
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